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GENERAL INFORMATION’ 


he i ‘BACKGROUND. 

“PLSIR was written by Dr R. oe Sandhu of the Ohio State, University. 

te: coding was based upen ‘the earlier “work _ Professor Wilson of the : 
"University of California at ; Berkley, but ‘modifications. were s made Ae | 
“improve t the program’ Ss: efficiency. 7 The © present version in use by AFEDL/- 


| FBC includes modifications by Mr . J “Muha to: account for orthotropic : 


meerteis and to permit bie: output, of “either maximum and minimum ‘stresses 


or strains, as -- 

| M2 2 PROGRAM DESCRIPTION 

-‘This program performs an elastic, plane etveedy finite element, 
structural analysis. It can handle linearly varying cherie) and pressure 


‘loads. “At che present » time it can handle qua aes and triangular 


elenents . up to three thousand erent ae ‘or grid points, a to gyeive mat- 


erials, up ‘to > eight temperature points for computing material properties. 
‘of each rear ak pressure acting on 1 up to three hundred elements. 
; The maximum semi-band width is fifty. | : 

The ‘analysis proceeds from 
vhere A is the stiffness matrix, x is the displacement matrix, and F is | 
he force matrix. ‘Fe | | 
3 Having determined the displacenents from Equation (1), the strains — 


an 


are assumed to follow from) - 


al. 


- where | B is the matrix linking the strains and the displacements, ard ¢ 
is the strain “patrix. 
Finally, taking the strains from Equation (2), ene: stresses are 


"found from 


‘where o is ee stress matrix and ( C is the material stiffness matrix. 
as... 4 tee fsotropic naterials, Equation 9) is. Hooke! Ss Law. 
ks 3. FUTURE WORK | | | as 
“As stated in Section i Ls. PLSTR has undergone a substantial modifi- | 
ae since being received by AFEDL/FBC. Future modifications now being 
 eaasaneee ae: using SPLINE interpolation pancei ue) 16 incorporate non- 
linear material properties, a plane strain option, the ability to tun 
| multiple load cases, and a failure criterion for developing margins of — 


safety. 


nes tip tag 7 ae: a . 7 ue : a nC) i" 


“PART i 


INPUT INSTRUCTIONS 


The Amput for PLSTR ee of sige logical cards. It must be 
“noted ‘that a “Logical card may consist of more than one physical card. 
Ose | For the remainder of this section, a Logical ee will be referred to 
simply as a aaa 
bee aa ‘The eight input cards needed to run PLSTR are f formatted as - £ollows: 


~ CARD 1: Title | 


“Columns Format Variable Explanation 
1-80. 8al0 |) «HED re a Any alphanumeric information | 
te ae | nr re. rn necessary to identify the problem 
CARD 2: Basic Information 7 a aa | 
Columns Format | Variable _ oe “Explation 
ae ae oe | _ NUMNP Se Number of grid points — 
6 - 10 ao a5 sa  NUMEL 7 ‘ | " Nunber of elements |. 
11 15. - 15 — NUMMAT na | | Number of materials 


16-20. 15 NuMPc) = -~CO™*~é<“«~*~*C«Nmbeerre atm : 
- ee See OD ee a i eae Gee Card 7) 


“OL - 30 F10.2 © ACELR' — : Meeelecation? in x-direction ee 
31 - 40° -F10.2 _ 3 ACELZ | | - Acceleration in y-direction 


“ge. abe 50 F102 Q - Reference (stress-free) ° 
| oe = ao | 3 pi a8. temperature | - 3 


CARD 3: Material Identification 
Colums Format Variable | wa _ Explanation 7 
ee 5 | eo MIYPE- i . Material Identification Number 


6 ~ 10 oe Loe a | NIC(MIYPE) | 7 Number of temperature cards: amy &,! 
Be BM RP ate ey eee we ode, ey — for material MIYPE (see Card 4). 


- Colums — 
“41-20. 


QA 25 - 


“CARD 3A: 


_ Columns 


1-10 — 


Format -- Variable eg ae 


‘F10.0 RO(MTYPE) 


15s NORTHO(MIYPE) 


| Material Orientation (Orthotr 
Format - Variable 


F10.0 Angle, 0 


opic Materials Onl 


= Explanation | 
Mass density for material MIYPE 


“Material type for material’ <6 .% 
MIYPE - | | a. 


0, ISOTROPIC MATERIAL = 


‘1, orthotropic material 


‘Explanation 


Orientation angle of material | 
‘longitudinal strength direction 
(fiber direction for dilamentary 
composites) in x-z plane 


_ Figure 2-1. Grid Geometry ~ 


ai 


GARD iA: Material Properties. (Isotropic Materials only) 
me Cone | ‘Format "Variable oe . 2 Explanation 
1-10) ‘Flo. 3 ORG, 1,MIYPE) | a: rT 
ic. Di 2302 F10.3° £(1,2,MIYPE) . Young's Modulus ae Coup erative: 1. 
21-30. ©¥F10.3 "BCL, 3)MIYPE) -_ Saint s Ratio at. aaa 
 % 31-40 oe. E(1,4,MIYPE) Poy of thermal expan~ 
oe. cat WA ae ee oe oe ee Mees - ' gion at ecmperature L 


} - _ NOTE: ” ihepeutr card: 4A £6t poche beuneteeare desired for material MTYPE, 
e | | a Css mepeee: NUMEG (HEZPE) times. 7 | | ; 


CARD 4B 4B: Witerial Pro eae (Orthotrop: ic Materials Onl 
| Columis’ Format ea Variable “ — Explanation | 


4-10. F10.0 THAT (MTYPE, 1) <_< Temperature 1 
41 O00. 81000" 7 £11 (MTYPE, 1) a : Lenigitudinal Youle! S “modulus at 
io as | ee | 7 | temperature | ner fe eS 
21-30  FIO.0° E22(MIYPE,I) | . ‘Transverse Young' Ss modulus ee. | 
t Fe, he a ; ge SS pPeReeraE I | a 


31 - 40. 10.0 <a G12(MTYPE, 1) ae Shear modulus in the 1, 2 plane 
_ | at temperature i ; 
“a = 50 F10.0 | -AMU12(MTYPE,1) Poisson’ Ss ceeicte oie ; 


temperature i 


51 - 60 | F10.0 _ AL(MTYPE,I) _ | - Tanet tad ined eoettictent ae 


: a » os a 4, aa > ae es wa Beem i expansion at oe a if 
| — 61-70. ©F10.0 | A2(MTYPE, I) ‘Transverse coefficient of thermal 
aa fo: g ee Ss =. oO expansion at Pompe tee tre. I. : 


, SO 71 = 80 - -F10.0 - (A12(MTYPE, I) an Shearing - soetricieae of ‘emai! 
ee tate Se Se es ae a oe | Fe expansion, 199 at a tats ae I 


a NOTE FOR (£12 (MTYPE, I) - At this time 12 has not been incorporated into poke 


Teal ATE Teh MN oe 2X tere xo 93 2 


the enalysis}, it has been included in the input with ‘possible future inclu- 


_sion in mind. For “orthotropic materials ae = ea 


ee at 


NOTES Repeac Card 4B for each t emperature. desired for material MIYPE, ~~ 


| oi. seas repeat Nunc (rr ees 


Ms ae; 


Cards” 3 aad. a are repeated NOMMAT times. 


~~ CARD 5. 5: crid Point Data : 


Essie | Format see ‘Variable ce a Explanation a 
Ps ae 15 : oN a ae - Grid point number as 


6-10 £F5.1 . ©. CODE(N) te eoandaee:. cgudicion rare on grid 
ee ee pudides)  .. i * oi poi N: 0, URN) | and. UZ(N) are 
- x and y loads’ |. = 
1, UR(N) is: x-displacenent and 
-UZ(N) is y-load . | | 
> 2, UR(N) is x- -load ana | uZQN) is 
_ y-displacement' 
3, UR(N) and UZ(N) are x and y 
| displacements | ion 


“A720, Oe RW) i x. “X-coordinate of grid point N 


“31 =930 F104 - 28) ae he aes | Y-coordinate of. erie point N 


ae 31 -.40 : F1O.4 UR(N) _ “X-load « or aicpieenene of orid 


point N(See CODE(N) above) 


41 ~ 50 F10.4 - -UZ(N) _ nl 44s load or displacement of “ha 
7. bg <P Cag - ae _ point N (See ene papOve) 


oe 60 7 F10.4- TON) 2 ad Temperature of erid point No 
Card 5 is repeated to ener bi teeeie ees system, ‘The repetition is 
performed, either for every grid point. or for those grid — required 
by the mesh generator within PLSTR cece Appendix on 7 | 
ea 6: ie etauea. | | 
= eowannis cere | - Variable - as Explanation 
1-5 im 1 ae ae er a - Element number . 


| 6-10 15 Per 8 IX(M,1) 00 o. First grid point for element M 


Wei, us * 5 IX(M,2)==~=——~—.—~«:s«SSe@cond grid point for element M. * 


16 - 20 Sy : Cn iX(M,3) U0 - ‘Third grid point for element M 
. 2125 I5 -IX(M, 4) ee Fourth Berd point for elenent M 


26 -- 30 Bee oe TXMAS) a Material identification. for Mg ae 
| | | “F* 7 J glenunn M , : 5 


as. for Card 5, Card ie. is ‘repeated either, for every element | or - for ‘those 


elements required by’ ‘PLSTR! s mesh generator Gee Appendix a). The 


. "sequencing of grid points | for an element is counter clockwise as shown _ 


"dm Figure” 2-2. 


Figure 2-2 “Element Grid Point Sequencing | 


_. (A) Quadrilateral oe @) Tetangle 
ioe ie ae ee ee A 


“NOTE: Four grid points must be input for each element types 

ie is: Bena Data a NUMPC is greater than zero) | 

‘Gelanae Format Variable Pa oo Explanation : 
4 : 5 « as: | _ IBC(L) tee pressure grid point 
ee 10 oe 7 gee) 2 Second pressure grid point: | 
a 20 F10.3 PR(L, 1) | =e Pressure. acting at TBC 


~— 


a - 30 - -F10. 3 PR(L, 2) ie Pressure acting at IBC. s 


| Card 7 is ‘repeated NUMPC times (See card 2). “The grid points are sequenced — 


“in a counter clockwise manner as you proceed! around the perimeter. of the ae 


fF 


“grid (See Figure : 2- 3). ‘A linear distribution is assumed along the ele- 


ment face. - 


Bch) 


 PR(L,1) 


Columns — 


Figure 2~3 Pressure Card Nomenclature 


——~ sB0(L) 


PR(L,2) 


; Output Selection 


Format  § Variable © a io Explanation — 


I5 NPRT “ee - Stress-Strain output flag 
ars is rn v 0, print for each element | 
El. No., x,y, x-stress, y-stress,- 


a an oS  yestress, max-stress, gun S Ete res 


and angle 


1, print for each element El. Now, 
X,Y, x-Stress, y-stress, xy-stress, . 
‘x-strain, y-strain, and xy-strain 3 


+ 


"APPENDIX AL 


"-PLSTR MESH GENERATOR a 


“The mesh 1 generator in | PLSTR is loacted between labeled statement 60 


“(ine 34) and labeled statement 190 (Line 68) ‘of the main program (See 


: appendix: B). ‘This generation Schone is useful in reducing the required 
ee: A number of physical puPae A cards when many handle. similarly sized ele~ 
een occur in the a, | | 
| “The “impact oF on ari cards needed is that, 1 if the N + si consecutive. 
sees points, M to M + N, are evenly spaced the variable, “CODE, “for each 
0 these points is zero, and 1 no loads are applied only ‘the erid points is 
: and M + N need be entered. The generator will 1 calculate — 
7 ie, es . *MEN-K,, . 
ta, 4 ia ag 
and 
an a a TYHN-Y,, 
, a 
a, fa we Then ‘ie mesh generator will assign 
eS ey Hae Gane 7 


ye 7 Tea Yea “Ta 4 


7 Temperature will be handled in the | same 2 manner, and ¢ CODE for alk generated : 


owes ethane im it Im, Set ee el Aad ASE en i ST ty 


points: will be set to zero. 


i Rated oa 


: Ce a ee oe an LP Paar eg ery YT MEE 
eee ea ama SLA sakiil rE canbe ON MERE eee EO been CAA 
. ? ‘ &, es : ; ' . 


eis 


The impact of the mesh generator on ‘the number ‘of element cards needed _ 


is similar to the impact on 1 grid cards. Tf the X a PS consecutive | 2 
“nents, M to M + x, “have erid points progressing byt ones, and the material 
ie the same for all N +. L elements, only the ‘elements M. cand M+ N need ibe 


Melee: The mesh oes ca assign 


IX(J-1,K) #1 J = MEL,MHs Ke1,40 7 


“and ee | 
'. Ix(J,5) = 18045) J MHL, MEN-1 
- ‘The following examples will illustrate the use and effect of the mesh ° 
generator. St ai 


: EXAMPLE 1. Take a a grid consisting of 64 grid points and ! 49° elements (See 7 
Figure AD). ‘In. this problen, assume material : 1 exists ‘from Y ai 0. 0 to 

G. 32, and material 2 exists from Y= 0. 32 to cp a Also assume that no 

"points are cou e cuarp ery 1 e@ os CODE is zero everywhere. Without: a mesh 


generator, 64 grid point cards and 49 element cards would be needed. 


The grid plate which must ae input for: this problem are ve 6; 8, 9, 14, 


- 16, 7, 22,24, 25, 30, 32; 3338, 40, 41,46, 48, 49, 54, 56, 57, 62, 64. ‘Thus, only 24 


| grid Hommre: need be input instead of 64. The elements car must be input 


are ‘4, 4, 5, ?, 8,11, 12, 14,15, 18, 19, Zig 22, 25, 26,28, 29, an, 33,35; 36, 39, 40, 42, 43, 


46, 47, 49. “Therefore, only 28 “elenents ‘need be input iapeenal of. 49. In : 


‘this problem, the mesh generator causes 61° fewer cards to’ be required. 


210. 


T 


[ev manola | 


Ll 


EXAWLE 2 2. “Take the sane gtid as in example 1. Assume, this time, ‘that 


bay. loads are > applied at grid points a to 8, , a zero x-displacenent, is | imposed 


ae grid points 57 to 64, and “only: one material is used. 


- The grid points which must be. input for this problen « are 18, 9; 14, 16 aA 7 
_- 22,24,25,30, 32, 33, 38,40,41, 46, 48,4954, 56,57- 6h. The elements which must | | 
be repel are a 7,8,14, 15, 21, 22, 28, 29. 35, 36, 42, 43, and 49. Since only 34 | 
| grid points and 14 elements ‘are input, 65 fewer sate cards are required 
“because of. the nesh g generator. | | eas | 
‘EXAMPLE 3 3. ‘Take a grid consisting of 100. grid points. end 1 162 2 elenents 


(See Figure A-2). - Aseame one material, ‘no loads, and no constraints. 


fc \ Without the mesh generator, 262 data cards would be’ required. 


The grid points which must be input are 1, ‘lo, do 205 21, ous 31, 40, ‘1, 50, 
ee ,60, 61, 70, 71, 80, 81, 90, 91, and 100. The elements which must be input 

_ are il, 9, 10, 18, 19, 27. 28, 305375 45, 46, 54, 55,63, 64, 725133845 82, 90; 91, 99, 100, 
108, 109, 47, 118, 126, 127313355130; 144, 145,153, 154, and 162. “Thus, only 56 : 


cards are required, instead of res 


990 
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eee in the x-direction. . A constant pressure of 100 psi is applied 
to dis owes surface. ee ee ee ree 7" 70°F, 
while the stress-free temperature of the structure is o°F, The gece 
core consists of two materials: as follows: | | 
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